Discrete wave turbulence of rotational capillary water waves 
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We study the discrete wave turbulent regime of capillary water waves with constant non-zero 
vorticity. The explicit Hamiltonian formulation and the corresponding coupling coefficient are ob- 
tained. We also present the construction and investigation of resonance clustering. Some physical 
implications of the obtained results are discussed. 



I. INTRODUCTION 

Capillary water waves are free surface water waves, 
the range of wavelengths being A < 17nim. In spite of 
their small amplitude, capillary waves play an important 
role in the dynamics of the ocean surface. Wind gen- 
erated capillary waves provide roughness on the ocean 
surface which then yields subsequently the a ppe arance of 
capillary-gravity and gravity waves [s^. [btI. |43| . On the 
other hand, capillary waves can be pumped by gravity 
waves [36', 'mI and therefore they also transport energy 
from bigger to smaller scales. 

The role of capillary water waves in the general theory 
of wave turbulence is invaluable. One can say that most 
new groundbreaking results and effects in (weak) wave 
turbulent systems have been first discovered for capillary 
water waves. Seminal work of Zakharov and Filonenko 
[SSj laid in 1967 the very foundation for what later be- 
came known as statistical (or kinetic) wave turbulence 
theory: the first power energy spectrum has been ob- 
tained for capillary water waves, as an exact solution of 
the wave kinetic equation. The wave kinetic equation 
describes the wave field evolution at the finite inertial in- 
terval < (fci, fc2) < 00 in fc-space. At this interval the 
energy of a wavesystem is supposed to be constant while 
forcing (usually in the small scales < kjor < fci) and 
dissipation (in the scales k2 < kdis) are balanced |5g |. 

Wave kinetic equations are deduced for the continu- 
ous A:-spectrum in Fourier space while in bounded or 
periodical systems the wave spectrum is discrete and 
described by resonance clustering |'29'|. An important 
fact first shown by Kartashova in [.25] is that indepen- 
dent resonance clusters might exist in the infinite spectral 
space and are nonlocal in fc-space (2^. [25|. Moreover, it 
was shown numerically (for atmospheric planetary waves) 
that increasing the wave amplitudes till some limiting 
magnitude yields no energy spreading from a cluster to 
all the modes of the spectral domain but rather goes 
along some "channels" formed by quasi-resonances, that 
is, resonances with small enough frequency broadening 
(Fig. 2, d^). Non-resonant modes with comparatively 
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big frequency broadening do not change their energy at 
the appropriate time scales (Fig. 3, [25|). 

A similar effect for irrotational capillary waves has 
been discovered by Pushkarev and Zakharov [s^ [s^ who 
investigated how the discreteness of the spectrum of ir- 
rotational capillary water waves influences the statistical 
properties of a weakly nonlinear wave field described by 
the kinetic equation. It was shown in a series of numer- 
ical simulations with the dynamical equations of motion 
(with energy pumping into small wavenumbers) that the 
wave spectrum consists of ^ 10^ excited harmonics which 
do not generate an energy cascade toward high wavenum- 
bers — frozen or fiuxless modes. "There is virtually 
no energy absorption associated with high- wavenumbers 
damping in this case" ([s^, p. 107). The quasi-resonance 
channels of the energy flux have also been observed and 
it was shown that the size of (quasi)-resonance clusters 
grows with increasing of the resonance broadening ((39j. 
p. 113). Two mechanisms for increasing the number of 
flux-generating modes have been identified: 1) growing 
the modes' amplitudes, and 2) enlarging the number of 
modes. 

The difference between the results of [2^ and [1^ 
is as follows. Irrotational capillary waves do not pos- 
sess exact resonances [2^ which yields the appearance of 
only fluxless modes; on the other hand, resonances do 
exist among atmospheric planetary waves. Accordingly, 
in [2^ two types of modes had been singled out — frozen 
and flux-generating. The specifics of the discrete wave 
turbulent regime is due to the fact that its time evolu- 
tion is governed by the dynamical equations, not by the 
kinetic equation [29] ; accordingly, phases of single modes 
are coherent and random phase approximation (necessary 
for deducing the wave kinetic equation) does not hold. 

The pioneering result of Zakharov et al. (6l| was 
the discovery of the mesoscopic wave turbulence regime 
which is characterized by the coexistence of discrete and 
kinetic regimes in numerical simulations modeling tur- 
bulence of gravity waves on the surface of a deep ideal 
incompressible fiuid. The mathematical justification of 
this co-existence was given in in the frame of the 
model of laminated wave turbulence. 

A peculiar property of irrotational capillary waves is 
that they do not have exact resonances in periodic boxes 
^^at, telai^r®tt)ilen(aieMhdifiSH|ete regime is exhibited in the 
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form of frozen turbulence, with distinct modes just keep- 
ing their initial energy. As it was shown quite recently, 
rotational capillary waves do have exact resonances, for 
specific positive magnitudes of the vorticity [8], which 
means that their discrete regime will have a different dy- 
namics. The simplest physical system where non-zero 
vorticity is important is that of tidal flows, which can be 
realistically modeled as two-dimensional flows with con- 
stant non-zero vorticity, with the sign of the vorticity 
distinguishing between ebb/tide [l6j . 

The main goal of our paper is to study the discrete 
wave turbulent regime of rotational capillary water waves 
with constant nonzero vorticity. In Section |lT] we give a 
brief overview of different turbulent regimes for three- 
wave resonance systems and we present known results 
on irrotational capillary waves (with zero vorticity) . De- 
tailed presentation of discrete wave turbulent regime is 
given in Section ITTll In Section HV] we provide some back- 
ground material on water waves with vorticity. Reso- 
nance clustering is constructed and studied in Section IVl 
A brief discussion concludes the paper. 

II. WAVE TURBULENT REGIMES 

To introduce the three generic wave turbulent regimes 
we need to distinguish between the notions of quasi- 
resonance and an approximate (non-resonant) interac- 
tion. Let us call a parameter S, 

|i±2±...±s| =<5() >0, (1) 

resonance width. Regarding the dispersion function j 
as a function of integer variables which are indices of 
Fourier harmonics (corresponding to a finite-size box) , it 
was shown by f28^ that the resonance width can not be ar- 
bitrary small. This means that if 6{) ^ 0, then 6{) > TZ{) 
where TZ{) is some positive number. Accordingly, ex- 
act resonances, quasi-resonances and approximate inter- 
actions can be introduced as follows: 

(1) exact resonances are solutions of ([T|) with 

60 = 0; (2) 

(2) quasi-resonances are solutions of ([!]) with 

o<d{)< no- (3) 

(3) approximate interactions are solutions of ([T]) with 

7^() < (5() < 7^ 

max 0, (4) 

where TZrnaxO defines the distance to the nearest reso- 
nance. The number TZO can easily be found numerically, 
for a given dispersion function, as in (46| . 

A. The model of laminated wave turbulence 

According to ©-([I]), discrete and continuous layers in 
a wave turbulent system are defined as follows. 



The discrete layer is formed by exact resonances and 
quasi-resonances governed by ([2]) and (jSj , while the con- 
tinuous layer is formed by approximate interactions given 
by Q. Four main distinctions between the time evolu- 
tion of these two layers have to be accounted for: 

1. dependence or independence on the initial condi- 
tions, 

2. locality or non-locality of interactions, 

3. existence or nonexistence of the inertial interval, 

4. random or coherent phases. 

The continuous layer 

This layer is described by the wave kinetic equations and 
its stationary solutions in the form of power energy spec- 
tra. The main necessary assumptions in kinetic wave 
turbulence theory are random phase approximation, lo- 
cality of interactions and existence of a finite inertial in- 
terval In the kinetic regime the time evolution of the 
wavefield does not depend on the details of the initial 
conditions and not even on the magnitude of the initial 
energy. Only the fact that the complete energy of the 
wave system is conserved in the inertial interval is im- 
portant. Interactions are local, that is, only wave vectors 
with lengths of order k are allowed to interact. More- 
over, interactions take place only at some finite interval, 
and no information is given about time evolution of the 
modes with wavenumbers lying outside this interval. 

The discrete layer 

This layer is governed by the set of dynamical systems 
describing resonance clusters; its time evolution depends 
crucially on the initial conditions. Modes' phases are 
coherent; the linear combination of the phases corre- 
sponding to the choice of resonance conditions (called 
dynamical phase of a resonant cluster) describes phases' 
time evolution (see Section UlI B[) . Resonant interactions 
are not local, that is, waves with substantially different 
fc-scales may form resonances. For instance, gravity- 
capillary waves with wave lengths of order of k and of 
the order of k^ can form a resonance; more examples are 
given in [3(| . Also the existence or absence of inertial in- 
terval is not important: resonances can occur all over the 
fc-spectrum. Moreover, in some wave systems waves with 
arbitrary wavelengths can form a resonance, for instance 
rotational capillary waves \§\ . We will discuss this aspect 
in detail in Section IVl 

Main aspects 

The main features of the model of laminated turbu- 
lence can be summarized as follows. 

• Independent clusters of resonantly interacting modes 
can exist all over k-spectrum. 

• The time evolution of the modes belonging to a clus- 
ter are not described by KZ-spectra, thus leaving the 



"gaps" in the spectra; their energetic behavior is de- 
scribed by finite (sometimes quite big) dynamical sys- 
tems. Depending on their form and/or the initial con- 
ditions, the energy flux within a cluster can be regular 
(integrable) or chaotic. 

• The places of the "gaps" in the spectral space are 
defined by the integer solutions of resonance conditions 
(e.g. see below). They are completely determined 
by the geometry of the wave system. 

• The size of the "gaps" depends on a) the form of dis- 
persion function, b) the number of modes forming mini- 
mal possible resonance, and c) sometimes — though not 
always — also on the size of the spectral domain. 

• The discrete layer can be observable all over fc-space 
and is characterized by the energy exchange among a 
finite number of resonant modes. The continuous layer 
can be observable within the inertial interval (fci, fc2) and 
is characterized by the energy transport over the fc-scales. 

Accordingly, three distinct wave turbulent regimes can 
be identified with substantially different time evolution, 
as shown schematically in Fig[T] Energetic behavior of 
the modes in the discrete regime is governed by a set of 
dynamical systems (see next section). 

Energy transport in the kinetic regime is described by 
the kinetic energy cascade. 
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FIG. 1: Color online. Schematical representation of three 
wave turbulent regimes. In the upper panel two co-existing 
layers of turbulence are shown; (black) solid arrows are di- 
rected to the wave turbulent regimes in which a layer can be 
observable (shown in the middle panel). (Red) dashed arrows 
going from the lower panel to the middle one show the evo- 
lutionary characteristics of the wavefield which manifest the 
corresponding regime. 



Ekm - fc-", (5) 

where > is a positive number depending on the dis- 
persion function and the number of interacting waves. 

Much less is known about the mesoscopic regime where 
both types of the wavefield time evolution can be de- 
tected simultaneously: energy exchange within resonance 
clusters and energy cascades. The existence of energy- 
cascading clusters has been first shown in fs^. The sim- 
plest dynamic cascade (without dissipation) has the form 

Edyn-p-'', (6) 

where p is a constant, < p < 1, n = 1,2,...,A^ and 
N is the number of cascading modes within a cluster. 
Dissipation can also be included, i.e. dynamic energy 
cascade is not affected by the existence or non-existence 
of the inertial interval. 

Whether an energy cascade in the mesoscopic regime 
has kinetic or dynamic origin can easily be checked in 
experimental data. 

In various physical systems either continuous or dis- 
crete or both layer(s) can be observable. For instance, in 
the laboratory experiments reported in fl7( only discrete 
dynamics has been identified. Coexistence of both types 
of the wavefield evolution has been detected in laboratory 
experiments, ^5l\, and in numerical simulations, ^6^]. In- 
clusion of additional physical parameters could yield the 
transition from the kinetic to the discrete or mesoscopic 
regime, |8j]. 



B. Parameter distinguishing among wave turbulent 
regimes 

In order to observe the various wave turbulent regimes 
in laboratory experiments and in numerical simulations 
one has to carefully make a few choices. 

First of all, it is necessary to define the range of al- 
lowed amplitude magnitude so that the waves belong to 
the "corridor" of weak nonlinearity. This means that 
the amplitudes should be big enough to leave the linear 
regime but should not be too big (to prevent strong tur- 
bulence). For water waves the wave steepness e, which 
is the characteristic ratio of the wave amplitude to the 
wavelength, is taken as a suitable small parameter, the 
usually choice being e = 0.1. 

Secondly, the choice of the spectral domain, i.e. the 
number of Fourier modes, is of the utmost importance. 
It is usually regarded as the main parameter allowing 
to distinguish between different wave turbulent regimes. 
Namely, the ratio of the characteristic wavelength A to 
the size of the experimental tank L is, as a rule, cho- 
sen as the appropriate parameter for estimating whether 
or not the waves "notice" the boundaries. For capillary 
water waves, if L/A > 50, the wave system can be re- 
garded as infinite and the dispersion function, both in 
circular and in rectangular tanks, can be taken in the 
usual form as fc'^/^ (M. Shats, private communication). 
The different choices explain why various laboratory ex- 
periments with capillary water waves demonstrate dif- 
ferent turbulent regimes. For instance, in [44.] discrete 
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modes have been identified while in [2^ the development 
of full broadband spectra is observed. 

Characteristic wavelengths 

It is tempting to use the parameter 

L/X (7) 

as a characteristic allowing to distinguish between dif- 
ferent wave turbulent regimes: taking a sufficiently large 
number of modes one will get the classical kinetic regime. 
However, the situation is more complicated. As it was 
shown in (57| . Fig. 4, low- frequency excitation (beginning 
with 25 Hz) generates distinct peaks of wave frequen- 
cies while increasing the excitation to 300 Hz yields an 
isotropic turbulent regime with a power law distribution 
of energies. Further increasing the excitation frequency 
generates an unexplained peak at 400 Hz. Moreover, cap- 
illary waves in Helium do notice the geometry of the ex- 
perimental facilities for L/ A ~ 30 — 300, the correspond- 
ing dispersion in circular tanks being described by Bessel 
functions (G. Kolmakov, private communication). In this 
case there exist nonlinear resonances and their manifes- 
tation has been recently established in laboratory exper- 
iments with capillary waves on the surface of superfluid 
Helium, [l| . The local maximum of the wave amplitudes 
is detected at frequencies of the order of the viscous cut- 
off frequency in the case when the surface is driven by a 
periodical low-frequency force. This means that the en- 
ergy is concentrated in a few discrete modes at the very 
end of the inertial interval (fci,fc2), in the beginning of 
the dissipative range. 

These novel results are paramount both from theoret- 
ical and practical point of view. This is the first exper- 
imental confirmation of the predictions of the model of 
laminated wave turbulence that resonant modes can ex- 
ist within and outside the inertial interval in the large 
fc-scales. 

The origin of these peaks (similar to rogue waves in 
the oceans) in the wave spectrum can be seen from the 
Fig|21 whether or not a rogue wave is generated depends 
both on the driving frequency drive and on the resonant 
frequencies in a wave system. Indeed, suppose the modes 
Ml, M2, M3 have frequencies 1,2,3 and form a resonance 
1-1-2=3. If drive ~ 3 (the modc with high frequency 
which is unstable), then the appearance of peaks with 
frequencies i and 2 is to be expected. If the modes 
M2, M^, M4 form a resonance with 2 — 3 + 4, then an 
excitation with drive ~ 2 yields peaks at 3 and 4. The 
importance of the high-frequency mode is discussed in 
the next section. 

Notice that rogue waves can appear within and out- 
side the inertial interval while discrete wave turbulent 
regimes are not connected with the existence of the in- 
ertial interval. The situation is quite general and does 
not depend on the type of waves under consideration. 
This is in agreement with the experimental observations 




FIG. 2: Color online. The energy power spectrum Ek is de- 
picted as the (grey) solid curve. Energies of resonant modes 
are shown by (red) vertical solid T-shaped lines having a 
string-like part. The energy of fluxless modes is shown by 
(blue) vertical sohd T-shaped lines. The empty circles in 
corresponding places of the energy power spectrum represent 
"gaps", indicating that for these modes the energy is given not 
by the solid curve but by the T-shaped lines. The horizontal 
and vertical axes show wavenumbers and energy correspond- 
ingly, while fci and ^2 denote the beginning and the end of 
the inertial interval. 

in superfluid Helium where rogue waves have been de- 
tected and the fact that the energy balance is nonlocal in 
/c-space was established [T9| . 

Characteristic resonance broadening 

The usual way to characterize the physical meaning of 
the resonance width 6 is to regard it as a shift in the 
resonant wave frequencies, cf. (40l. l6l| . yielding reso- 
nance broadening in a wide wave spectrum with a large 
number of modes. 

To characterize the broadening in a three- wave system 
governed by ([T^ one can choose as an appropriate pa- 
rameter the inverse nonlinear oscillation time of resonant 
modes 

r^' « \ZBk\. (8) 

In discrete regimes the phases of individual modes are co- 
herent; however if the width | A^Sg | of the high-frequency 
mode in a resonant triad becomes substantially larger 
that TjT^, the coherence is lost — this being a necessary 
condition for the occurrence of the kinetic regime. One 
concludes then that 

• The discrete regime corresponds to 

|A,B3l«T,7\ (9) 

• The kinetic regime corresponds to 

|A,B3l»T,7\ (10) 

• The mesoscopic regime corresponds to 

|A,B3l^T,7^ (11) 
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In the case when exact resonances are absent and only 
approximate interactions have to be accounted for, the 
estimates (|9l)- (|TT|) can be made general (cumulative) by 
the substitution \Ak\ instead of jA sgj, where |Ak| is de- 
fined as a broadening of a mode with wave vector k. The 
fact that the generation of the kinetic regime occurs via 
spectral broadening of discrete harmonics has been es- 
tablished in ^40j for irrotational capillary waves. 

There are three main reasons for the absence of exact 
resonances in a wave turbulent system with decay type of 
dispersion function: a. the resonance conditions do not 
have any solutions ( ^ fc^/^); b. the resonance conditions 
do not have solutions of a specific form (for example, for 
^ solutions exist in a square but not in a majority 
of rectangular domains); c. all resonances are formed in 
such a way that at least one mode in each primary clus- 
ter lies outside the inertial interval. In this case other 
modes will become "frozen" ([H,!!^) and no resonance 
occurs in the inertial range of the wavenumbers. In the 
case when exact resonances are not absent, one has to 
be very careful while introducing any cumulative esti- 
mate: for approximate interactions the upper boundary 
T^maxi), which depends on k, has to be included. The 
problem of introducing a suitable cumulative parameter 
for distinguishing among various wave turbulent regimes 
is contingent on future extensive research. 

III. DISCRETE REGIME 

The discrete wave turbulent regime is characterized by 
resonance clusters which can be formed by a big number 
of connected triads. The dynamical system correspond- 
ing to a cluster formed by N triads can be written out 
explicitly by coupling the N systems for a triad 

Biy = ZjBl\^B^j,B2\j = ZjB^jB^^j, 

B^Ij ^ ~ZjBi\jB2\j, (12) 

with J = a, 6, c, ... and equaling the appropriate B^\j. For 
instance, the dynamical system of the first two-triad clus- 
ter shown in Fig (3] reads 

{Bl\a = ■^aS2|a^3|a, B2\a = ZaB^^B^\a^ 
B^a = —ZaBi\aB2\a, -Si|f, = ZbB*|j^i33|f,, 
B2\b = ZbBl^^^B^\i, iJgih = -Zbi?i|hB2|6, ^ 
Bl\a = Bi\f, , i?3|a = B2\b 

{Bl\a = ZaB*^^B3ia + ^fc-B3|a-^3|b, 
B2\a = ZaBl^^B^la, ^^g^ 
B3\a = -ZaBi\aB2\a + ZhB^^^B^^^, 
B-i\b = -ZbBi\aBz\a, 

The system (|13p is unique, up to the change of indices 
1 O 2. On the other hand, the form of the system 
changes depending on the fact whether or not the con- 
necting mode is a high-frequency mode in one or both or 
in no triads. 



The high-frequency modes B^a and B^, are 
called active, or A-modes, and the other modes 
{Bi\^, B2\a, Bi\jy, B2\b) are called passive, or P-modes (s^ . 
because due to the Hasselman's criterion of instability 
[2^ ^-modes are unstable under infinitesimal excitation 
and P-modes are neutrally stable. The cluster dynamics 
is then defined by three possible connection types within 
the cluster, namely AA~, AP~ and PP— connections. In 
the system ([T^ the connection Pi|a = Pi|t is of PP- type 
and P3|<j = P21& is of AP-type. 

The graphical representation of an arbitrary resonance 
cluster in the form of a A'TJ-diagram {NR for nonlinear 
resonance), first introduced in [29.], allows us to keep this 
dynamical information. The vertices of a 7Vi?-diagram 
are triangles corresponding to resonant triads and the 
half-edges drawn as bold and dashed lines denote A- and 
P-modes correspondingly. Examples of A^i?-diagrams are 
shown in Fig. [3] 

Each mode of a cluster generates a gap in the KZ- 
spectrum. The maximal possible number of gaps can 
also be seen from the form of the A^P-diagram, though 
the actual number can be smaller, in the case when two 
different wavevectors fci and k2 have the same length. 




FIG. 3: Color online. Af7?-diagrams for clusters generating 
maximally 3, 4, 5, 6 and 7 gaps (from left to right and from 
up to down). 

An illustrative example can be found in [4] where 
laboratory experiments with two-dimensional gravity- 
capillary waves are described. By analyzing the mode 
frequencies in the measured data, only five different fre- 
quencies were identified: 10, 15, 25, 35 and 60 Hz. How- 
ever, theoretical consideration allowed to conclude that 
in fact seven distinct modes take part in the nonlinear 
resonant interactions, and the corresponding resonance 
cluster has the form of a chain formed by three connected 
triads. The amplitudes and frequencies were identified as 

Ai o 60Hz, A2 ^ 35Hz, ^3 o 25Hz, 
Ai ^ 25Hz, A^ o lOHz, Aq ^ 25Hz, Ar o 15Hz, 

and resonance conditions for frequencies as 

1=2 + 3, 2 = 4 + 5j 6 = 5 + 7 (14) 

"The waves A3, A^ and Aq all have the same frequency 
(25 Hz), but they must have different wavevectors, in or- 
der to satisfy the kinematic resonance conditions. We 
assume that the mechanical means to generate a test 
wave at 25 Hz also generates perturbative waves in other 
directions at 25 Hz" ([i|, p.70). The A^P-diagram of this 
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FIG. 4: Color online. Af7?-diagram of the resonance cluster ap- 
pearing in the laboratory experiments with gravity-capillary 
water waves reported in [4^. 



cluster is shown in FigHJ the cluster is a 3-chain with one 
PA- and one PP-connections. The dynamical system cor- 
responding to the cluster has been solved numerically: 
for all calculations the dynamical phase of the initially 
excited triad was set to tt, while during the experiments 
the phases of the modes have not been measured (D. 
Henderson, private communication). This yielded quali- 
tative agreement of the results of numerical simulations 
with measured data but higher magnitudes of observed 
amplitudes. 

The evolution of the amplitudes and phases in the dis- 
crete turbulent regime is discussed in Section IIII Al and 
in Section IIII Bl correspondingly. 



A. Evolution of the amplitudes 

In Fig[2]we show the difference in dynamics of discrete 
modes of capillary waves with and without vorticity. The 
energies of fluxless modes in the frozen turbulence regime 
for irrotational capillary waves are shown as (blue) bold 
vertical lines and the energies of resonant modes appear- 
ing due to the non-zero constant vorticity are depicted by 
(red) bold lines with a spring part inside. The "springs" 
show symbolically that the magnitudes of energies can 
change in time, similarly to the process of compressing 
or pulling an elastic string. Notice that this situation 
is quite general and does not depend on the wave type 
(see [25.] for numerical simulations of atmospheric plan- 
etary waves — in Fig.l the periodical time evolution of 
resonant modes is shown and in Fig. 3 the fluxless modes 
demonstrate no time changes in the magnitudes of am- 
plitudes). 

The analogy with a string is quite helpful for the un- 
derstanding of the dynamics of resonant modes. A string 
can be compressed or pulled (without destruction) only 
in a finite range of the applied forces; similarly, the range 
of amplitudes' changing is also finite. Indeed, the sim- 
plest way to show it is just to look at the smallest possible 
resonant cluster — a triad. The dynamical system for a 
resonant triad (fT^ has the Hamiltonian 



Ht - lm{B,B2B;) , 



(15) 



and two Manley-Rowe constants of motion [35(| in the 
form 



(16) 



providing the integrability of IT^ in terms of the Jaco- 
bian elliptic functions sn, cn and dn [55^. Notice that 
(jl2p has 6 real variables (real and imaginary parts of the 



amplitudes Bj, j = 1, 2, 3) and only 3 conservation laws 
P^ . (fTB|) which, generally speaking, are not enough for 
integrability. However, if we rewrite it in the standard 
amplitude-phase representation Bj — Cj exp{i6j) as 



Ci ~ ZC2C3 cos(p, 

C2 = ZC1C3 cos if, 
6*3 = —ZC1C2 COS (p 



(17) 



C2 ^ C's 



it becomes clear immediately that in fact we only have 4 
independent variables: three real amplitudes Cj and one 
dynamical phase 



ip = Ul 



(18) 



according to the choice of resonance conditions. In this 
section we will concentrate on the wave amplitudes while 
the importance of dynamical phases will be discussed in 
the Section IIIIBl Making use of the addition theorems 
sn^-|-cn^ = 1 and dn+/i'^sn = 1 (here /i is a known func- 
tion of conservation laws, = ii{H, I23, 113) and of the 
amplitude-phase representation Bj = Cj exp(i6j), one 
can easily get explicit expressions for the squared ampli- 
tudes: 



Cj{t) = aij + a2j ■ dn'^^Q;3j -t-Z, a^j 



(19) 



with the coefficients aij, a^j being explicit functions of 
-^13, I23 and Ht, for j = 1, 2, 3, [lH. The energy of each 
mode Ej is proportional to the square of its amplitude, 
Ej(t) ^ Cj(t) ^ dn^, i.e. its changing range is finite and 
defined by the initial conditions as < dn^(t) < 1. 

There is no analytical expression for the amplitudes 
(energies) of a cluster formed by two or more triads. On 
the contrary, as it was shown in 0, H^, HO], the corre- 
sponding dynamical systems are integrable only in a few 
exceptional cases. For instance, clusters formed by N 
triads all having one common mode and with all con- 
nections of AA- or PP-type are integrable for arbitrary 
initial conditions if ZijZj = 1/2,1 or 2. Some clusters 
are known to be integrable for arbitrary coupling coeffi- 
cients but only for some specific initial conditions. On 
the other hand, already the smallest possible cluster of 
two connected triads can present chaotic behavior (30j . 
though the modes energies are still bounded as the cor- 
responding dynamical systems are energy conserving. 

In the case of a generic cluster, the time evolution of 
the amplitudes has to be studied numerically and de- 
pends crucially on the initial conditions. In order to de- 
crease the number of degrees of freedom one has to use 
the fact that each triad has 3 integrals of motion given 
by P^(|16p . i.e. isolated triads have 2A^ conservation 
laws and A^ Hamiltonians. A cluster formed by A^ triads 
with n connections has (2A'^ — n) conservation laws and 
one Hamiltonian. The integrability of a generic cluster 
depends on the connection types within a cluster and on 
the ratios of coupling coefficients Zj. Some results of 
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numerical simulations with two-triad clusters (fixed con- 
nection type, various coupling coefficients) can be found 
in 0. 

B. Evolution of the dynamical phases 

As one can see from P7)) . the time evolution of a res- 
onant triad does not depend on the individual phases of 
the resonant modes but on the dynamical phase cor- 
responding to the resonance conditions. In this case, the 
closed expression for the evolution of dynamical phase, 
first found in [3l|, reads 

ip{t) = sign((y5o)arccot(^fci J'l^fca {t-to),fi^^ , (20) 

where 

and 

J^{x, ^) = sn{x, ^) cn(a:, ^) dn(a;, ^), 

while (^0 is the initial dynamical phase and the constants 
IJ,,K{p),T are known explicitly as functions of /13, I23 
and Ht- 

The time evolution of a generic cluster again does not 
depend on the individual phases but on the dynamical 
phases; the corresponding equations for the dynamical 
phases have a form similar to the last equation of ([T7| . 
For instance, in case of two triads a and b connected via 




FIG. 5: Af7?-diagram for a two-triad cluster with one PP- 
connection. 

one PP-connection (the A^i?-diagram is shown in Fig (5]), 
the equations for the dynamical phases 

fa = 6l\a + ()2\a — ()3\a, <^h = + ^'2|6 — ^3|fc (21) 

read 

^a = ~H{C^' + C-^l-C-p, (22) 
= -i?(Cr' + C-^l - C-^l). (23) 

where the joint mode is denoted by Ci and H is the 
Hamiltonian 

H = Ci {ZaC2\aC3\a Sm (fla + ZbC2\bC3\b Sin fb) ■ (24) 

The effect of the dynamical phases on the evolution of 
the amplitudes has been studied numerically in It 
turned out that if initially ipin = for a triad, then it 
stays zero and the changing range of amplitudes is max- 
imal. If ifin is very small but non-zero, say, (fin = 0.01, 
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FIG. 6: Upper panel. Color online. Plots of the dynamical 
phases as functions of time, for each frame: (pa{t) is (red) 
solid and ipbit) is (black) dashed, with (pa,in ~ fb,in = 7r/2. 
Lower panel. Plots of the squared amplitudes, as functions 
of time. The connecting mode Ci is blue (the upper one), 
C2\a and Csia are green (in the middle), while C2\b and Csij, 
are purple (at the bottom). To facilitate the view, C||i,,C3ji, 
and Ci are shifted upwards. The horizontal axis stands for 
non-dimensional time and the vertical axis for the phase in 
the upper panel and for the squared amplitudes in the lower 
panel. 



then this range is substantially diminished: depending 
on the initial energy distribution within a triad it can de- 
crease 10 times and more, being minimal for ipin = tt/2. 
For the cluster shown in FigIS] the maximal range of the 
amplitudes is observed for ipa,in = fb,in = 0. 

An interesting phenomenon has been uncovered in the 
case when at least one of the phases ipa and tpb is initially 
non-zero: a new time-scale is clearly observable where the 
amplitudes of the resonant modes are modulated by the 
evolution of the dynamical phases. The characteristic 
form of evolution is shown in FiglHlfor fa — ft — t^/^, 
i- 

This is a manifestation of one more crucial distinction 
between discrete and statistical wave turbulent regimes: 
while in the first regime phases are locked, in the latter 
they are supposed to be random. The effect of phase 
randomization among three- wave interactions of capillary 
waves has been studied recently experimentally in [40j . 

The transition from a coherent-phase to a random- 
phase system occurs above some excitation threshold 
when coherent wave harmonics broaden spectrally. The 
waves were excited parametrically in the range of the 
shaker frequencies 40 3500 Hz and the modulation in- 
stability of capillary waves, which can be approximated 
by the squared secant function sech^, has been estab- 
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lished. Experiments show that the gradual development 
of the wave continuum occurs due to the spectral broad- 
ening. Similar effects have been previously observed in 
liquid Helium [JTl and in spin waves in ferrites [ssj. 

However, the form of the resulting KZ-spectrum does 
not allow to conclude automatically whether the energy 
cascade is due to three- or four- wave interactions, and 
both explanations can be found in the literature. A di- 
rect way to clarify the matter would be to estimate the 
periods of the excited modes for the case of three- and 
four-wave interactions at the corresponding time scale. 
The experimental confirmation can then be obtained by 
measuring the corresponding wave fields |42|]. 



IV. FLOWS WITH CONSTANT VORTICITY 

A rotational water wave is a wave in which the un- 
derlying fluid flow exhibits non-zero vorticity. Physi- 
cally rotational waves describe wave-current interactions: 
a uniform current is described by zero vorticity (irrota- 
tional flow) 0, [HI , while a linearly sheared current has 
constant non-zero vorticity [18]. The irrotational flow 
setting is appropriate for waves generated by a distant 
storm and entering a region of water in uniform flow 
[s^liwhile tidal flows are modelled by constant vortic- 
ity |16l |. The existence of periodic gravity wave trains 
with small and large amplitudes (T2| . was established, 
and qualitative properties of such waves were studied 
d, M, M, M,M,M,M,\M- In the case of capiUary- 
or gravity-capillary waves the rigorous theory is limited 
to the small-amplitude regime [52], although some par- 
tial results pertaining to the existence of large-amplitude 
waves have recently been obtained [53]. Note that al- 
ready the small-amplitude theory is more involved in the 
presence of surface tension, including such phenomena as 
Wilton ripples [56| in which two harmonics with the ratio 
1 : 2 are in resonance. 

An important new result on rotational capillary waves 
is the theorem on the dimension of flows with constant 
vorticity [sj: 

Theorem. Capillary wave trains can propagate at the 
free surface of a layer of water with a flat bed in a flow 
of constant non-zero vorticity only if the flow is two- 
dimensional. 

This theorem allows us to consider for flows of constant 
vorticity only two-dimensional flows propagating in the 
x-direction. For irrotational flows the existence of a ve- 
locity potential enables the transformation of the govern- 
ing equations for capillary water waves to a Hamiltonian 
system expressed solely in terms of the free surface and 
of the restriction of the velocity potential to the free sur- 
face The absence of a velocity potential for non-zero 
vorticities complicates the analysis considerably and one 
can not expect results of this type. 

However, for flows of constant vorticity the remarkable 
preservation of the main features of the above theory (for 



irrotational flows) can be established by introducing a 
generalized velocity potential. In particular there is a 
Hamiltonian formulation in terms of two scalar variables, 
one of which is the free surface elevation [sj] . In contrast 
to the three-dimensional waves, it is convenient to regard 
the y-axis as the vertical axis, with the waves propagating 
in the a;-direction. The problem can be written in terms 
of a generalized velocity potential, satisfying 



A(j) = 0, — oo < y < J], 
|V0| 0, y^-oo. 



(25) 



with (px — u -\- and (j)y — v, |51j] . In this coordinate 
system the free surface oscillates about the mean water 
level y = 0. 

Introducing the new variable C, with 



n 1 



(26) 



where f — 4>\y=r] is the restriction of the generalized ve- 
locity potential to the free surface, allows to obtain a 
canonical Hamiltonian system in the form |51|j: 



. _5H ■ SH 



(27) 



where the notation 5F/5u is used for the variational 
derivative of a functional F with respect to the variable 
u. 

Passing to Fourier variables can be seen as a change of 
variables and this transforms Hamilton's equations into 



. _ SH ■ SH 

'54 ^Vk 



(28) 



for k e Z*. 

The change of variables 



with 



(29) 



linearizes the quadratic Hamiltonian and transforms 
Hamilton's equations into 



In the new variables we have that 

H2=J2 ^ik)\ak\ 
feez* 

where 



(30) 



(31) 



oj{k) = -^sgnk-^^<j\k\^-^^, (32) 
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is the dispersion relation for two-dimensional 27r-periodic 
capillary waves with constant vorticity [52] . 

The non resonant terms in the Hamiltonian can be 
eliminated up to any desired order by the normal form 
transformation, i.e. we can replace the Hamiltonian func- 
tion H = H2 + H3 + Rihy H = H2 + iJf " + ^4. Here 
H^'^^ is the resonant part of H^, 



E 



S{ki, k2, k'i){ak^ak2a*k.i + c.c), 



k-^—k\-\-k2 
uj{k3)—uj{k\)+uj{k2) 



(33) 

with S{ki, k2, fcs) being known functions of wavenumbers 
kj first found in f^. 

Let us now restrict our attention to a resonant triad 
with resonance conditions taken in the form 



i:j(fc3) = a;(fci) -I- ^(^2), ^3 = ^1+^2 



(34) 



If we approximate the Hamiltonian using the quadratic 
and cubic terms, Hamiltonian equations take the form 



di + zwifli — —iSiki, ^2,^3)0203, 
d2 + iuj2a2 = —iS{ki, fc2, ^3)0103, 
0.3 + iuJzaz = ~iS{ki, fc2, ^3)0102. 



Setting aj 
the standard form 



, these equations transform into 



Bi = ZB2*S3, B2 = ZBIB3, B3 = -ZB1B2, (35) 
where the coupling coefficient reads, [9j, 

n. 



kik2k3 

TTW 10)2 0)3 



UJ1UJ2 



-UJ3 



(36) 



Note that in this formula the variables kj , luj , dij and 
n are not independent, since they are related by the res- 
onance conditions ([M)) . All possible magnitudes of the 
constant non-zero vorticity generating resonances, i.e. so- 
lutions of ([M)) . can be directly computed, [8], as 



nik,,k2) = 

kik2 i9kl + 9/c| + 14/cifc2)(fci + k2)-^/^ 
6 ' y/{6kf + 15k'fk2 + 22k'ik'i + Ibkik'i + 6fc|) 



(37) 



for arbitrary fci, ^2, k^ satisfying (j34p . 

It follows from p7p the magnitude of a positive vortic- 
ity triggering a resonance can not be too small. The min- 
imal magnitude of the constant non-zero vorticity which 
generates resonances is 



(38) 



The main conclusion one can deduce from (|37p is that 
any two rotational capillary waves with arbitrary wave- 
lengths can form a resonance for a suitable magnitude 



of vorticity. In other words, any vorticity computed by 
([57)1 will generate an isolated triad; its dynamics has been 
discussed in detail in Section IIIII 

However, in laboratory experiments the magnitude of 
vorticity can only be controlled within some non-zero 
£ error corresponding to the available accuracy of the 
experiments. For capillary waves e ~ 10~^ is easy to 
achieve and more refitments allow to reach e ~ lO""*^ in 
some cases. In Section |V] we construct the resonance 
clustering of rotational capillary waves for various ap- 
proximate magnitudes of the vorticity. 



V. RESONANCE CLUSTERING 

To see how the resonance clustering depends on the 
available accuracy e, let us define the e-vicinity of vortic- 
ity as 



fir. 



(39) 



and let us construct examples of clustering for a few dif- 
ferent choices of e. 

Notice that since the vorticity depends linearly on ^/a^ 
the previously defined e does not depend on a and conse- 
quently the resonance clustering constructed below will 
be the same for capillary waves in an arbitrary liquid. 
However, even for the same liquid, say water, the magni- 
tude of (T depends on the temperature and atmospheric 
pressure while performing experiments. All numerical 
simulations discussed below have been performed in the 
spectral domain fci,fc2 < 100. 

• 10-8 < e < 10-5. 

We began our numerical computations with e = 10^^ 
which does not generate any other clusters but one iso- 
lated triad for each magnitude of the approximate vor- 
ticity |r2±£|. Increasing e from 10^* to 10^^ still leaves 
us with isolated triads. 

• 10-5 < £ < 10-"*. 

For this values of e, four new clusters appear, each 
formed by two connected triads: 



'(20, 94, 114), (24, 70, 94); 

(17, 71, 88), (15, 88, 103); 

(11, 83, 94), (12, 71, 83); 
[(10, 47, 57), (12, 35, 47). 



(40) 



Each line in (|40p consists of two resonance triads and the 
exact magnitudes of the generating vorticities are given 
by ([371), say 17(20,94) and f2(24, 70). These magnitudes 
do not coincide, 17(20, 94) 7^ 17(24, 70), however 



0(20,94) - 0(24,70) 
' 0(20,94) ' 



< e. 



(41) 



All four two-triad clusters have the same NR-diagram 
shown in Fig 17] No results on the integr ability of cor- 
responding dynamical systems are presently known, and 
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namical system reads 



FIG. 7: NR-diagram for clusters appearing for e < 10 



we expect chaotic energy exchange among the modes of 
such a cluster. 

• 10-4 <e< 10-3. 

The structure of resonance clustering is substantially 
richer in this case. For instance, 83 different vortici- 
ties can generate two-cluster clusters, among those most 
clusters have AP-connections but also clusters with AA- 
connections appear, for instance 



(50, 50, 100), (49, 51, 100); 
(47, 48, 95), (46, 49, 95); 
(44, 44, 88), (43, 45, 88) 



(42) 



and others. 



FIG. 8: NR-diagram for clusters described by (|42|) . 

The corresponding NR-diagram is shown in FiglHl The 
explicit form of the dynamical system for a cluster con- 
sisting of a triad a and a triad b connected via two high- 
frequency modes, i?3|a — Sai&j reads 



Bl\a 


= ^aB^\^B3\a , 


Bl\b = ZbB*^i^B3\a , 




< B2\a 


= ZaBl^^B^^a , 


B2\b = ZbBl^i^B^^a , 


(43) 


.■Bala 


= -ZaBi\aB2\a 


— ZbBi\bB2\b- 





The integrability of ([321) has been studied in 0, with 
the following results. The system ([^5)) is integrable for 
arbitrary initial conditions^ if Za/Zb is equal to 1, 2 or 
1/2, [i^. In the case of arbitrary Za/Zb, the integrability 
can be proven for some specific initial conditions [sj . 

Beside two-triad clusters we also observe the appear- 
ance of various structures formed by three, four, five and 
more connected triads; some examples are given below. 

The cluster of three triads 

(79, 80, 159), (78, 81, 159), (77, 82, 159) (44) 
has the NR-diagram shown in FiglHl accordingly, its dy- 




FIG. 9: NR-diagram for the clusters given by (|44p (on the 
left) and by (gSJ (on the right). 



Bl\a 


^ ^a^2\a^5\a ' 


B2\a = ZaB*^^B3\a. : 


Bl\b 


= ZbB'^^fjB-i\a J 


B2\b = ZbBl^f^B^ia , 


Bi\c 


= ZcB2\i.B^a , 


B2\c = ZcB^^^B^^a , 


B3.\a 


= -ZaBi\aB2\a 


— ZbBi\bB2\b — ZcBi\cB2\c 



(45) 



where the indices a, b and c are taken for three triads; 
all connections within this cluster are of AA-type, i.e. 



B 



3|a 



B 



3|b 



B 



3|c 



The cluster of four triads 

(48, 48, 96), (47, 49, 96), (28, 96, 124), (46, 50, 96) (46) 

has three AA-connections and one AP-connection. The 
NR-diagram (see FiglHl) ^-'i^ the dynamical system can 
be easily constructed but no results on the integrability 
of the corresponding dynamical system are known. 
• 10-3 <£< 10-^ 

The size of the clusters grows exponentially with the 
growth of e so that for e — some vorticities generate 
clusters formed by a few hundreds to few thousands of 
connected triads. The maximal cluster in the studied 
spectral domain consists of about 4000 triads with more 
than 33.000 connections among them. 

It would be plausible to assume that in this situation 
no regular patterns will be observable for generic initial 
conditions. On the other hand, a special choice of initial 
conditions — excitation of a P-mode — might produce a 
fluxless (frozen) regime. However, one has to construct 
the resonance clustering and to check whether or not 
the "parasite" frequency of electronic equipment used 
during experiments generates a resonance (see [4J for 
more detail). 

Last but not least. Let us rewrite (|57|l as 

^{ki,k2)=(T^'^ ■^{ki,k2) (47) 

where 

kik2 i9kf + 9fc| + 14fcifc2)(fci + k2)-^^^ 



^{kiM 



yj&i&kl + 15fc3fc2 + 22fc2fc2 4- l5A:ifc3 -|- 6fc|) 



does not depend on the properties of fluid: these are 
absorbed in the coefhcient a. The expression for Cj{k) 
can also be rewritten as 



(48) 



with (I A; 1 3 -I- Cl^ /4) not depending on a. Correspondingly, 
the interaction coefficient Z has the form 



Z = . z 



(49) 



with Z depending only on the wave numbers ki,k2,k^. 
This means that the cluster integrability defined by the 
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ratios of the corresponding coupling coefficients [4^ [s^l 
also does not depend on the properties of fluid a while 

Z,JZy=~Z,rZy (50) 

where Zj and Zy do not depend on a and the indexes j, j' 
correspond to the j-th and the j'-th triad in a cluster. 
However, the magnitude of the vorticity $1 generating 
the corresponding cluster will depend, of course, on the 
properties of a fluid as a function of a. 

Indeed, a is ratio of the coefficient of surface tension 
tr to the liquid density p, i.e. a = oj and is different 
even for the same liquid if experiments are performed 
under different conditions (such as a change of the ambi- 
ent temperature). For instance, for water with p = 998 
kg/m^ and standard pressure, a — 72.14 mN/m if the 
temperature of water is 25° C and a = 75.09 mN/m if 
the temperature is 5 C. Accordingly, cr — 7.23 • 10^^ and 
a = 7.52 • 10"^ for these two cases. 



VI. DISCUSSION 

Our main conclusions can be formulated as follows: 

• Among the rotational capillary waves with constant 
non-zero vorticity resonances can occur only if 1) the flow 
is two-dimensional, and 2) the magnitude of vorticity is 
larger than ilmin given by ([55]) . 

• Two arbitrary rotational one-dimensional capillary 
waves can form a resonance only for appropriate magni- 
tudes of the constant non-zero vorticity, with the exact 
magnitude given by ([37]). Thus a chosen magnitude of 
the vorticity generates one isolated resonance triad with 
a periodic energy exchange among the modes of the triad 
and the magnitudes of the resulting amplitudes depend 
on the initial dynamical phase (|20p . In a laboratory ex- 
periment where initially the A-mode of the triad is ex- 
cited we expect the appearance of some regular patterns 
on the surface of the liquid. 

• If some non-zero experimental accuracy is taken 
into account, clusters of more complicated structure oc- 
cur. Their time evolution can be regular (integrable) or 
chaotic, depending on the ratios of coupling coefficients 
and sometimes also on the initial conditions. For an ar- 
bitrary cluster its dynamical system can be written out 
explicitly and solved numerically, thus predicting the re- 
sults of a laboratory experiment beforehand. 

• The coupling coefficient for the corresponding dy- 
namical system is given explicitly by p6|) . Note that 



from of ([5^ it follows that in fact this expression gives 
the correct coupling coefficient for rotational capillary 
waves, not only in water but in arbitrary liquid. 

• Last but not least. The discrete turbulent dynamics 
shown in Figl2] occurs for very general types of flow mo- 
tions and is also valid for wave systems possessing 4- wave 
resonances though the corresponding theoretical study is 
more involved. Due to Hasselman's criterion for 4-wave 
systems [l^l, there is no analog to the A-mode in this 
case, while the excitation of any one mode in a generic 
resonant quartet does not yield an energy flow within a 
quartet. At least two resonant modes have to be excited. 
The only exception is given by the quartets of the form 
1+2 = 23, which can be regarded essentially as a 3- 
wave resonance and all the results presented above can 
be applied directly, with the A-mode having frequency 
23. The dynamical system for a generic quartet with 
frequency resonance condition 

a;(ki)+a;(k2) =a;(k3)+a;(k4), (51) 

is also integrable in terms of Jacobian elliptic functions 
(45| . The time evolution of a quartet is defined by the 
fact whether or not two initially excited modes belong 
to the same side of the (fSTj) or to the different sides. 
Correspondingly, the pairs (1,2) and (3,4) are called 1- 
pairs, while the pairs (1, 3), (1, 4), (2, 3), (2, 4) are called 2- 
pairs. The excitation of a 2-pair again does not generate 
an energy exchange among the resonant modes, while 
the energy pumping into a 1-pair will produce it under 
some conditions (see [lOl, Ch. 4, for more details). The 
A^i?-diagrams keep track of this dynamical information 
in the form of the edges [29] . Applying these results for 
gravity water waves (work in progress) we can see that 
in this case the generation of rogue waves demands more 
conditions to be fulfilled than in a 3-wave system. 
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